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Abstract 

An equation, based on the free volume of a liquid solvent, was derived 
via dimensional analysis, to predict binary diffusion coefficients. The 
equation assumed that interaction between the solute and liquid solvent 
molecules followed a Lennard-Jones [6-12] potential. The equation was 
compared to other diffusivity equations and was found to give good results 
over the temperature range examined. 

Introduction 

The prediction of binary liquid diffusivities has a long and varied 
history. As early as 1905, the Stokes-Einstein equation was derived from a 
hydrodynamic model (1)— a large solute sphere diffusing in a solvent 
continuum. 


d ab =kt '/ 

For definition of symbols, see the listing of notations. 
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In 1955, Wilke and Chang (2), also using a hydrodynamic model, derived 


the following semi -empirical expression: 


= 7.4xlO- 8 U^ B ) 1/2 T/(|i B v2 s ) 


( 2 ) 


where, X B is an association number of the solvent: 2.6 for water, 1.9 for 

methanol, 1.5 for ethanol, and 1.0 for unassociated solvents. 

Starting in 1941, Eyring (3) and his co-workers (4) derived the 
following equation using absolute rate theory: 


D„ = [ (KT) / (E a v. b ) ] (N/V b ) exp (E^-E^/RT) 


( 3 ) 


This equation is based on a quasi -crystal line liquid structure model where the 
diffusing molecules "jump" through layers of solvent molecules. The geometric 
parameter, E A , equals the number of nearest molecules, in a plane, around a 
central molecule. In using this equation, Eyring and co-workers assumed that 
E A = 6 and that E^ B = b Dm so that equation 3 reduced to: 


= {KT/6\i g ) (n/v b ) ^ 


( 4 ) 
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Akgerman and Gainer (5) in 1972, building on Eyring's work, proposed the 
following equation: 


D M = (KT/ E a \i b ) WV b ) 1 / 3 (MW b /MW a ) 1/z exp(£^ - E D jRT> 


( 5 ) 


In this equation, the geometric parameter is given by: 


e a = 6 {v A /v B ) l '« 


( 6 ) 


and, unlike Eyring, they did not assume that E^ B = E Dja> . 

Prediction of diffusion coefficients using equations based on a 
hydrodynamic model has been successful when the solute molecules are large and 
spherical, and the solvent mblecules considerably smaller in size. The Eyring 
equation, based on an absolute rate model, has been found to predict diffusion 
coefficients too high by a factor of six. The Akgerman-Gainer equation has 
been found to give better values than the Eyring equation. 

All the equations cited (Equations 1, 2, 3, and 5) utilize the solvent 
viscosity (a transport property) to calculate the diffusion coefficients 
(another transport property). Naturally, it was highly desirable to develop a 
model to calculate liquid diffusivities without using viscosity. In 1959, 
Cohen and Turnbull (6) derived the following relationship between the self 
diffusion coefficient in a model liquid of "hard spheres" and the "free 
volume" of the liquid: 
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(7) 


D = A 

They theorized that the liquid free volume (V F ) redistributed itself 
leading to the formation of "holes" sufficiently large (V*) for diffusion to 
occur. The use of equation 7 requires the estimation of V F , V*, and the 
constants A and 7 which depend on the geometry and size of the liquid 
molecules. 

The concept of relating liquid free volume to diffusivity is an 
interesting one, but modeling the interaction of the liquid molecules using a 
hard sphere potential (which incorporates only a repulsive force) apparently 
leads to the necessity of using a "hole" volume. The purpose of this work was 
to investigate the possibility of formulating a diffusivity equation based on 
the liquid free volume where the molecules interact according to the more 
realistic Lennard-Jones [6-12] potential. The use of the Lennard-Jones [6-12] 
potential (which incorporates both repulsive and attractive forces) should 
eliminate the need for the "hole" volume. An equation capable of predicting 
both self and binary diffusivities was sought. 

Formulation of Free Volume Equation 

The initial formulation of the diffusivity equation assumed that 
diffusion coefficients should depend on the temperature and free volume of the 
liquid solvent. 



D, 


AS 


f(T,V„) 


( 8 ) 
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The diffusion of solute molecules through solvent molecules also depends 
on the interaction between the molecules and the molecular size and shape of 
the molecules. In this case, the molecular interaction was modeled using the 
Lennard-Jones [6-12] potential which has two associated parameters --a 
distance parameter (cr) , and energy parameter (e/K) . For solute molecule A 
interacting with solvent molecule B, the following empirical combining laws 
were used (7): 


= i ( ° A 


+ Ofl) 


(9) 


(c/K)ab ■ [U/IOaU/iOb] 1 ^ 


( 10 ) 


One can then calculate Q which is a dimensionless function of 

d ao 

T/(e/K) AB ; q d is derived from the kinetic theory of Chapman and Enskog, and 
is a measure of the deviation of the Lennard-Jones molecules from rigid-sphere 
behavior: q =1 for rigid spherical molecules. Values of o, (e/K), and 

Q have been tabulated. The empirical combining laws give good estimates 

d m 

for nonpolar molecular pairs. 

The molecular size was taken into account by empirically averaging the 
molecular weights of the solute and solvent: 


Mff M = <.m a +m b ) (n) 

Including the Lennard-Jones parameters and the average molecular weight 
into the initial diffusivity relation yields the following: 
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^AB Tf VfB* ®A8* 


( 12 ) 


Dimensional analysis was then used to derive the following dimensionally 
correct equation: 



where G is a dimensionless function of Q^. At this point, G was determined 

empirically. This was done by correlating equation 13 with the experimental 
diffusivity value of argon diffusing in carbon tetrachloride at 298K. The 
argon-carbon tetrachloride system is considered to approach the ideal 
conditions of a lennard-Jones [6-12] model--nonpolar, nonhydrogen bonded 
spherical molecules. Letting G = Q'^ (the simplest form for G was sought) 

correlated the experimental result well. 

The functional form for G was derived for the temperature equal to 298K. 
For temperatures other than 298K the exponent of q should change. Rather 

than manipulate the exponent, equation 13 was scaled by multiplying it by 
(T/298) n . The value of n was found by correlating equation 13 to the 
experimental diffusivity value or argon diffusing in carbon tetrachloride at 
273K. Letting n=3 resulted in a diffusivity value that deviated by less than 
1% from the experimental value (5); this results in equation 13 having a T 3 * 
temperature dependence. 

The liquid free volume may be calculated approximately by one of two 
methods. Glasstone and co-workers (3) derived the following equation for the 
calculation of liquid free volume: 
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where a is a geometric packing factor for the molecules of the liquid. It is 
equal to 2 for regular arrays of hard spheres. They also derived the 
following equation for the calculation of free volume: 



where V is the volume of a molecule in the liquid given by: 


v = 


mw b 
(p b») 


( 16 ) 


U Liq and U vapor are the velocity of sound in the solvent liquid phase and in 
the solvent vapor phase, respectively. Velocity of sound data for several 
liquids was obtained from the CRC Handbook (8). Calculated free volumes, 
based on velocity for some liquids are: 

V F (benzene) = 4.07 x 10' 25 cm 3 at 302. 6K 

V F (carbon tetrachloride) = 4.43 x 10' 25 cm 3 at 298K 

V F (water) = 6.75 x 10‘ 25 cm 3 at 298K 

The following example illustrates the use of equation 13 for the 
calculation of the diffusion coefficient for benzene, A, in carbon 
tetrachloride, B: 

R = 8.314*10 7 (g-cm 2 )/(sec 2 -gmole-K) 
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MW a = 78 g/gmole MW B = 153.8 g/gmole 

<j a = 5.349*10' 8 cm a B = 5.947*10’ 8 cm 

(e/K) A = 412.3 K (e/K) B = 322.7 K 

T = 298 K 

V FB = 4.43M0' 25 cm 3 
MW ab = 1/2 (MW a + MW b ) = 116 g/gmole 
a AB = 1/2 (a A + a B ) = 5.648*10' 8 cm 
(e/K) AB - ( (e/K) A * (e/K) B ) 1/2 = 364.8 K 
T/(e/K) AB = 0.82 4 fl D(AB) = 1.59 

D ab = 1/ (59) 2 * ( (8 . 3 1 4*10 7 ) * (298) / ( 1 16) ) 1/2 * ( (4.43*10 -25 ) / (5 . 648* 10" 8 ) ) 1/2 
D ab = 1.62*10* 5 cm 2 /sec 

Equation 13 overpredicted the experimental value by 5% (9). 

Table 1 lists the experimentally determined binary diffusion 
coefficients for a number of substances dissolved in different organic 
liquids. The calculated values using equation 13 are also listed, as well as 
the % deviation from the experimental values (5, 9) (the quantities in 
parenthesis). The calculated values and % deviation are also given using the 
Stokes-Einstein (S.E.), the Wilke-Chang (W.C.), and the Eyring, and the 
Akgerman-Gainer (A.G.) equations. 

The calculated free volume results differ from the experimental results 
by an overall absolute average value of 13%. The absolute average deviations 
for the other equations are 71% (S.E.), 21% (W.C.), 61% (Eyring), and 47% 
(A.G.). Comparing only the free volume results, note that the diffusivity 
values for C0 2 dissolved in C 2 H 5 0H and for N 2 dissolved in C 6 H 6 have the 
greatest deviation from the experimental values (-58% and +27%, respectively.) 
These results are not surprising if one considers that equation 13 
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These results are not surprising if one considers that equation 13 
incorporated Lennard- Jones parameters that are derived from a model that 
considered the interacting molecules to be nonpolar spheres with no hydrogen 
bonding. 

This point was explored further by calculating the diffusion 
coefficients of several different solutes in water— a polar, hydrogen-bonding 
molecule. The calculations are tabulated in Table 2. Equation 13 predictions 
are very poor compared to the other four equations. The free-volume overall 
absolute average deviation from the experimental values is 173%, whereas the 
deviations for the other equations are 46% (S.E.), 9% (W.C.), 30% (Eyring), 
and 10% (A.G.). In order to improve on equation 13 predictions, it was 
decided to use a polar potential, the Stockmayer potential, for water. 

The Stockmayer potential has three associated parameters— a distance (a) 
and energy (e/K) parameters (which are not equal numerically to the Lennard- 
Jones parameters), and a dimensionless polarity parameter (t*) which is a 
function of the molecule's dipole moment. These three parameters have been 
tabulated for a number of polar molecules (10). 

In order to incorporate the polarity parameter t , equation 13 was 
modified to the following: 


Dab 


(1-fc*) ( 


298 



( 


RTVfb ) 


1/2 


(17) 


The following example illustrates the use of equation 17: 

Binary diffusion coefficient calculation for C0 2 (A) in water (B) 
MW a = 44 g/gmole MW B = 18 g/gmole 

a A = 3 . 94* 10 -8 cm ff B = 2.52*10' 8 cm 
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(e/K) B = 775 K 


(e/K) A = 195 K 
T = 298 K 

V FB = 6.75*10" 25 cm 3 
t B * - 0.7 
MW ab =31 g/gmole 
a AB = 3.23*10' 8 cm 
(e/K) AB = 389 K 
T/(e/K) AB = 0.77 4 fi D(AB) = 1.64 

D ab = [(l-0.7)/(l. 64) 2 ] * ( (8 .3 14*10 7 ) * (298) / (31) ) 1/2 * ( (6 . 75*10 25 ) / 

(3 . 23*10 -8 ) ) 1/2 
D AB = 1 .44*10' 5 cm 2 / sec 
Experimental value = 1.91*10' 5 cm 2 /sec (5) 

Equation 17 underpredicted the experimental value by 25%. 

Table 3 lists the free volume diffusion coefficients, for various 
solutes dissolved in water, using equation 17. Comparison with Table 2 shows 
that using the Stockmayer potential for water results in calculated 
diffusivity values (overall absolute average deviation of 28%) that are in 
better agreement with the reported experimental values. The diffusivity value 
was also calculated for C0 2 in ethanol (C 2 H 5 0H) using the Stockmayer 
parameters for ethanol. Comparison of Table 3 to Table 1 again shows that use 
of the Stockmayer parameters for ethanol results in a better prediction (-17% 
deviation compared to a -58% deviation.) 

Equation 17 was finally used to calculate diffusivity values for C0 2 in 
ethanol and in water over a small temperature range. These values are 
tabulated in Tables 4 and 5. The calculated diffusivity values were also 
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graphed against the experimental values and these results are shown in Figures 
1 and 2. 

For C0 2 in ethanol, equation 17 overpredicts diffusivity values from 10% 
to 17%, whereas the other four equations underpredict from 41% to 78%. 

For C0 2 in water, equation 17 predictions are better than the Stokes- 
Einstein predictions and comparable to the Eyring predictions. Both the 
Wilke-Chang and Akgerman-Gainer equations gave better predictions than the 
free volume equation. 

Conclusions 

A new semi -empirical equation was formulated, via dimensional analysis, 
to predict both self and binary liquid diffusivities. The equation is based 
on the free volume of the liquid solvent and on the Lennard-Jones [6-12] 
molecular interaction between the solute and solvent molecules. Unlike other 
widely used diffusivity equations, equation 13 does not incorporate the liquid 
solvent viscosity (a transport property), and unlike equation 7, it does not 
need the "hole" volume. Equation 13 is able to predict both self and binary 
diffusivities with good results (over the temperature range considered) when 
compared with four other widely used diffusivity equations. This work also 
suggests that a similar equation may be derived to predict the viscosity of 
liquids. 
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Notations 

D self diffusion coefficient 

D ab diffusion coefficient of solute A and solvent B 

E a number of B molecules around a central A molecule on the same plane 

E D (ab) activation energy for diffusion 

E D activation energy for viscosity of solvent B 

p B 

G dimensionless function of q. 

U AB 
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AH 

un vapor 

K 

mw a 

mw b 

N 


r 

R 

Ra 

T 


U Liq 


u 


vapor 


FB 

* 


X B 

a 

(c/K) 


a 


P B 


enthalpy of evaporation 
Boltzman constant 
molecular weight of solute A 
molecular weight of solvent B 
Avogadro's number 
distance 

universal gas constant 
radius of solute A 
absolute temperature 

velocity of sound in solvent liquid phase 

velocity of sound in solvent vapor phase 

molar volume of solute A 

molar volume of solvent B 

free volume of liquid B 

hole volume 

association number of solvent B 
geometric packing factor 
Lennard-Jones energy parameter 
Lennard-Jones parameter distance 
viscosity of solvent B 
density of solvent B 
Lennard-Jones dimensionless function 
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TABLE 1 



% Deviation = 


Exp 1 1 value - calc, value 
exp'l value 


This definition used in Reference 5. 
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TABLE 2 


Solute 

A 

Solvent 

B 

Temp. 

(K) 

Exp'l D ab 

x 10 5 
(cm 2 /sec) 

Eq. 13 
D ab x 10 5 

S.E. 

W.C. 

Eyri ng 

A. G* 

ci 2 

h 2 o 

298 

1.45 

| 

■PWIB 

iKHlW 

Wfim 

1.65 

(-14) 

2.44 

(-68) 

1.60 

(-10) 

so 2 

h 2 o 

298 

1.94 

2.98 

(-54) 

0.95 

(51) 

1.75 

(10) 

2.44 

(-26) 

1 

wBSM 

o 2 

h 2 o 

298 

2.35 

WBBM 

1 

1.15 

(51) 

2.44 

(-4) 

2.44 

(-4) 

2.42 

(-3) 

n 2 

h 2 o 

298 

2.01 

1 

wBMM 

1.07 

(47) 

2.17 

(-8) 

2.45 

(-22) 

2.41 

(-20) 

CM 

o 

o 

h 2 o 

298 

1.91 

4.61 

(-141) 

1.03 

(46) 

2.04 

(-7) 

2.45 

(-28) 

2.03 

(-6) 


*Lennard-Jones (6-12) potential used for water. 
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TABLE 3 


Solute 

A 

Solvent 

B 

Temp 

(K) 

Exp'l 0 AB 

x 10 5 
(cnr/sec) 

Eq. 17 
D ab x 10 5 

S.E. 

W.C. 

Eyring 

■ 

Cl 2 

h 2 o 

298 

1.45 


wrm 

1 

1.65 

(-14) 

2.44 

(-68) 

1.60 

(-10) 

so 2 

h 2 o 

298 

1.94 

1.63 

(16) 

■8M 

1 

mm 

2.44 

(-26) 

EEI 

wtiSm 

o 2 

h 2 o 

298 

2.35 

2.24 

(5) 

1.15 

(51) 

2.44 

(-4) 

2.44 

(-4) 

2.42 

(-3) 

n 2 

h 2 o 

298 

2.01 


1 

Kgfl 

2.17 

(-8) 

2.45 

(-22) 

2.41 

(-20) 

co 2 

h 2 o 

298 

1.91 


■cm 

2.04 

(-7) 

2.45 

(-28) 

2.03 

(-6) 

o 

O 

ro 

c 2 h 5 oh 

298 

3.42 

4.01 

(-17) 

0.86 

(75) 

EH 

1.37 

(60) 

1.81 

(47) 


*Stockmayer Potential used for water and for ethanol (C 2 H 5 0H) 
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TABLE 4 


Solute 

A 

Solvent 

B 

Temp. 

00 

Exp 1 1 D ab 
xlO 5 

(cm 2 /sec) 

Eq. 17 
D»„ x 10 5 

S.E. 

w.c. 

Eyring 

A . G * 

O 

o 

ro 

c 2 h 5 oh 

279.4 

2.45 

2.85 

(-16) 

1 

WESm 

1 

its 

f!f!l 

1 

1.23 

(50) 

co 2 

c 2 h 5 oh 

283 

2.78 

3.05 

(-10) 

0.61 

(78) 

1.47 

(47) 

II 

WmM 

■KB 

MBS 

C0 2 

c 2 h 5 oh 

288 

3.02 

3.36 

(-11) 

0.70 

(77) 

m 

1.09 

(64) 

1.48 

(51) 

co 2 

c 2 h 5 oh 

298 

3.42 

4.01 

(-17) 

0.86 

(75) 

2.05 

(40) 

1.37 

(60) 

1.81 

(47) 

co 2 

c 2 h 5 oh 

303 

3.84 

4.44 

(-16) 

1 

mBm 

2.27 

(41) 

| 

wmm 

2.04 

(47) 


*Stockmayer potential used for ethanol (C 2 H 5 0H) . 
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TABLE 5 


Solute 

A 

Solvent 

B 

Temp. 

(K) 

Exp'l D ab 

x 10 5 
(cnr/sec) 

Eq. 17 
^AB X I® 5 

S.E. 

W. C. 

Eyring 

m 

co 2 

h 2 o 

2.83 

1.25 

nm 

m 

1 

mem 

1.58 

(-26) 

1.34 

(-7) 

co 2 

h 2 o 

288 

1.40 

mm u 

HLm 

0.78 

(44) 

1.55 

(-11) 

1.85 

(-32) 

1.55 

(-11) 

co 2 

h 2 o 

293 

1.70 

I 

w 

1 

1.80 

(-6) 

2.13 

(-25) 

1.80 

(-6) 

co 2 

h 2 o 

298 

1.91 

Mgm 

1.03 

(46) 

2.04 

(-7) 

2.43 

(-27) 

2.03 

(-6) 

CM 

o 

o 

h 2 o 

303 

2.27 

1 

umm 

m 

2.32 

(-2) 

2.75 

(-21) 

2.29 

(-D 


h 2 o 

313 


1.84 

(33) 

iffll 

H 

3.46 

(-25) 

2.83 

(-2) 

■9 

h 2 o 

323 

3.24 

2.18 

(33) 

1.81 

(44) 

3.60 

(-11) 

4.28 

(-32) 

1 

Bn 

co 2 

h 2 o 

338 

4.29 

1 

msm 

2.40 

(44) 

4.72 

(-10) 

5.62 

(-31) 

EE9I 

IBM 


*Stockmayer Potential used for water. 
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Diffusivity vs Temperature 



Figure 2. Comparison of Experimental and Free Volume Binary Liquid 
Diffusion Coefficients for the CC^-Water System 
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